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Magd E.Kahil ^ ^ 3 
Abstract 

The problem of motion for different test particles , charged and spinning objects 
of constant spinning tensor in different versions of the bi-metric theory of gravity is 
obtained by deriving their corresponding path and path deviation equations, using 
a modified Bazanski in the presence of Riemannian geometry. This method enables 
us to find the path and path deviation equations of different objects orbiting very 
strong gravitational fields. 


1 Bi-metric Theories: A Brief Introduction 

General relativity is considered a landmark in history of science of being, during the last 
century, as a pivotal non-linear theory of gravity [1]. Yet, some problems have been re¬ 
mained unsolved due to applying Riemannian geometry in its explanation. One of these 
difficulties is related to the law of conservation of energy and momentum [2], which is dealt 
with considering the metrical tensor as flat one at great distances from the gravitational 
source. Rosen [3,4] introduced a remedy to this problem by proposing two different met¬ 
rics representing the gravitational source and giving a curved space and describing 
a physical one expressing an inertial frame and becoming a flat space. Using this assump¬ 
tion, it can be found that the held equations of Einstein imply to a theory of gravitation 
in hat space . A problem arises because the pseudo-tensor quantities in Orthodox General 
Relativity, then turns not to be preserve its tensor character [2]. This led Yalmoz [5] to 
examine a new class of solutions for the held equations of Rosen’s theory of gravitation to 
solve the difficulty of dealing with the hat metric. But, a slight problem has emerged in 
that, the speed of light is no longer constant as conhrmed experimentally in the realm of 
Special Relativity. Such a problem is counted to be a virtue using Mohat’s approach of a 
bi-metric theory of gravity for a variable speed of light, as it helps to reveal the puzzle of 
the dark energy problem. This approach may be expressed by means of two metrics are 
likened to each other in terms of gradients of scalar /biscalar helds to explain the rapid 
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expansion of galaxies is due to change of speed of light from one epoch to another that 
i.e. dark energy can be interpreted This type of theories are expressing how the inflation 
scenario of the universe is due to Bi-metric version of variable speed of light [7] . Some 
applications using the Moffat formalism of bi-metric theory of gravity are explaining the 
causal description of quantum entanglement [8] and testing the propagation of neutrinos 
using OPERA experiment [9]. 

Another problem in GR is inability to explain the rotation curves of spiral galaxies. 
Milgram [10] proposed a specihc treatment by performing a modified Newtonian Dynam¬ 
ics paradigm (MOND)that was able to explain its causality apart from appealing to the 
dark matter problem. A new step has been taken successively to extend MOND to be 
expressed in terms of bi-metric theory to become BIMOND [11] by having two held equa¬ 
tions describing matter and twin-matter[12] may be used for examining the existence of 
gravitational waves and explaining two interacting 4D membranes[13]. Recently, Hassan- 
Rosen [14] developed an extension of the present bi-metric theory using the concept of 
bi-gravity, which used two metrics describing their gravitational helds. In doing so, they 
discard the earlier bi-metric theories metrics in which one metric describes gravity while 
the other is physical. Accordingly, a new massive gravity theory for spin-2 , free from 
ghosts, has been obtained [15]. Several applications of these types of theories are viable in 
dealing with obtaining held equations for very strong gravitational helds of neutron stars 
[16], which suggests the possibility to examine black holes and super-massive black holes 
e.g Sgr A* by studying the stability of objects orbiting in these helds. This hypothesis is 
essentially to be examined. 

2 Motion in Bi-metric Theories 

Theories of gravity may help to explain the functioning of diherent objects. For example, 
Rosen[3] obtained the equation of motion bi-metric theory of gravitation for a test particle. 
These equations were solved by Isreaht[17] in order to examine their behavior in the case 
of small velocities and weak helds using a post Newtonian approximation . 

Also, Falik and Opher [18] used the bi-metric theory of gravity to hnd the held equa¬ 
tions associated with spinning neutron stars as an example of a strong gravitational held 
. This achievement opens the way to examine the motion of diherent charged objects and 
spinning ones in the presence of a strong gravitational held as dehned by Bi-metric type 
theories. This led us to obtain the corresponding path and path deviation of diherent 
objects , such as test particles, charged particle, spinning objects. These results are de¬ 
termined by introducing a Lagrangian with a specihc feature for obtaining the path and 
path deviations. Therefore, the key role to this approach is based to obtain path and path 
deviation equations for each of these objects using a specihed Lagrangian for each case. 
The aim of our study is obtaining equations of motion for objects using bi-metric theory 
originated from two metrics in one stands for gravitational held and the other dehnes 
physical matter or both represent gravity. 
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2.1 Path and Path Deviation Equations: The Bazanski Ap¬ 
proach 

Geodesic and geodesic deviation equations can be obtained from the following Bazanski 
Lagrangian [19]: 

L = ( 1 ) 

where is the metric tensor [/" is a unit tangent vector to the geodesic , its deviation 
vector, S is the parameter characterizing the geodesic and geodesic deviation ^ is a 
covariant derivative with respect to g^j_v If one takes the variation with respect to , 


we get 

dr 

(2) 

and 


(3) 


d dL dU^ 

• = gaa,pU UP + g^a • 

dS d'^^ dS 

(4) 


Substituting from the above relations into the Euler-Lagrange equation: 


d dL dL 
dS 


(5) 


and taking into account the following condition: 


dfj.i'-tp — 0 - 


( 6 ) 


We obtain after some manipulations: 

dU°' 

= 0, (7) 

where is the Levi-Civita affine connection. Also, following the same technique can be 
applied to obtain the variation with respect to obtain the geodesic deviation equations: 


dL 








( 8 ) 


and 


d f)T 

dL 




d-jf " ■ DS 

Substituting in its corresponding Euler-Lagrange equation: 


d dL dL 
dS dU^ dx^ 


(9) 

( 10 ) 


( 11 ) 
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i.e. 


+ {9^Kf6^„U“),U'’ - (9, 


TJ^^ _ 


+ = 0 


The resultant equation is not tonsorially covariant unless, if one substitutes the following 
quantity 


DU^ 

DS 


0 , 


in it, and provided that the Riemann curvature is dehned as 




■pQ; _ pQ; j_ pe pa pe pa 

^ ^ pup ^ ^ pp^ ev ^ pu^ ep 


then, after some manipulations we obtain: 

_ — /?« \T>7 7'//5p 

jjg2 - ^ • 


( 12 ) 


The above method has been applied in different geometries than the Riemannian one 
e.g. non-Riemannian geometries admitting non-vanishing curvature and torsion tensors 
simultaneously [20-22], This approach helps to implement the concept of geometrization 
to include not only physics but also biological epidemic curves [23] as well as economic 
complex systems in terms of information geometry [24], Moreover, the Bazanski La- 
grangian has been modihed to describe the path equation of charged object to take the 
following form [25]; 

r)\i>4 p 

L = (13) 

Db m 

If we take the variation with respect to following the similar steps as in (5) we obtain 


dU^ 

~dS 




—FW 

m ■ 


(14) 


where is an electromagnetic tensor, ^ the ratio between charge to mass of an object. 
And, on taking the variation with respect to 17“ , the equation is not tonsorially covariant 


till we impose the following conditions 
(^\DU^ = e_ 

DS ^ 

{‘^)F^u-,p + Fpp.-,u + Fpp^u = 0 . 

Consequently, after some manipulations, we obtain its corresponding deviation equation: 


DS^ 


e 

+ —(F“—— + F 
m ■ Ds 


'.pup'- 


v\p 




(15) 


Moreover, for non processing spinning objects, equations of spin and spin deviation are 
obtained by applying action principle on the following Lagrangian : 


nu//? 1 

Db 2m 


(16) 


where is a spin tensor of a spinning object. By taking variation with respect to the 
4/“ we obtain 


717“ 

~d^ 





a 

.pup 


SPUR 19 


(17) 
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And taking the variation with respect to t/" , and using the following condition con¬ 
ditions: 

(1) w = , 

( 2 ) — 0 

to obtain its deviation equation: 


^ + Rp>,sfpr<i‘- + (18) 

Also, for deriving path equations for spinning charged object [27] we take the variation 
with respect to 'h" on the following lagrangian 


L 


Dm 




m 






fiDpa 




(19) 


to obtain 


dU^ 

Us 


+ r 


fiiy 


m ■ 2m 'p P 


( 20 ) 


And taking the variation with respect to 77“ , providing the following consitions: 

(1) w = iFfUf + 

(2) -F)ti/;p + Dp^-i, -|- Fpp.y = 0. 

(3) + Rip^ + = 0 

we obtain its corresponding deviation equation 


e Dm^ 1 

- = jjpupmp A —IF""-h F" upmp)A -F" jjpjjpmp 


Moreover, for a spinning object with precession we modihed Bazanski Lagrangian [29] 

nca/3 r)\S)0 1 

L = 9o.p{mU- + + -Rap^sS^'U^m- (22) 

to obtain equation of a spinning object by taking the variation with respect to the devi¬ 
ation vector '!/“ 

n n Qa/3 1 

= -Rl^pSPPUP (23) 

And taking the variation with respect to 77“ , and using the following condition con¬ 
ditions: 

(1) W = hFAysS-'^U^, 

(2) Rp^g + R2p^ + R'^sis = 0 

to obtain its deviation equation: 


Ds^ 


n Q12/3 p c\p 

R%pUP{mUP + Up^W + g--g.x{mU^ + 


\(R‘ 


gyp ^ ^ 

ppp Ds 


+ R^p^^SP.^^pUPmP + 


(24) 
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However, the spin precession eqnation can not be obtained from the above Lagrangian 
bnt dne to some relations between and its nnit tangent vector P" 


pi^U^ = (mU>^ + Ug )U’' 

U o 

(25) 

and 

DS‘'^ 

P-U^^ = {rnU’^ + Up )UP 

J-J O 

(26) 


Thus, after some manipulations, regarding that = 6^ one obtains the equation of 

spin precession; 


= 2{pf^U‘' - P'^U^). (27) 


2.2 Path and Path Deviation Equations in Weyl Geometry 

It is well known that in Weyl geometry the gravitational potential tensor is associated 
with such a scalar held. From this perspective one can dehne a combined gravitation 
potential tensor in the following manner [30]: 

^ (28) 

Where g^i, is the Weyl gravitational potential, and 0 a scalar held, which may give raise to 
introdnce disformal transformation of any gravitational theory having two metrics dehned 
in the following way [31] 

= [Ag^^ + (29) 

where A and A are arbitrary constants . 


Thus, in this type of geometry it can be dehned its corresponding affine connection to 
become : 

+ -g°'\ga54',g + g8g(t),a - gga4>,&) (30) 

In order to obtain the geodesic eqnation , we introdnce the following Lagrangian: 


L = (31) 

If we take the variation with respect to \h“, we obtain the path equation: 

rUja 

= 0, (32) 

And its corresponding deviation equation is determined by applying the commntation 
relation between the two parameters in a similar way as explained in (2.1) is obtained by 
taking the variation with respect to IP^ to become: 


DS^ 


R 


tiup 




(33) 
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where 


_ ■po; _ pa _l_ per pa _ per pa 

(ii) Dixon-like Equation for spinning charged objects of Weyl geometry: 

Similarly , we can obtain the Dixon-like path equation as dehned in Weyl geometry by 
suggesting the following Lagrangian 

df/" - e 1 - 

(34) 

and its corresponding deviation equation becomes: 


- e 1 - 


DS^ 


m 


DS 


^ + R%xS";:U-''t>’ + Ef;,y,S’'^U“'S^) (35) 


2m 


2.3 Path and Path Deviation Equations of MOND 

In this part, it is worth mentioning the path and path deviation of modihed Newtonian 
Dynamics (MOND) paradigm due to its vital role in explaining the vague regions due 
to dark matter problem that are unknown by Newtonian/Einsteinain formulations i.e. 
revealing the nature of rotation curves of spiral galaxies. Accordingly, some authors have 
studied motion of a test in MOND [32], This has led us to apply the Bazanski method 
in order to obtain the path and path deviation equations for any test particle related to 
this paradigm by suggesting the following Lagrangian 

1 

L = (36) 

DS m 

where 0 is a Newtonian potential added due to MOND . 

Applying the same approach of the Bazanski method in (2.1) by taking the variation with 
respect to T" we obtain its path equations 

+ = (37) 

Also, equations of geodesic deviation are obtained by taking the variation with respect to 
f/" on the Lagrangian (36) with taking into consideration that 


DU^ 

DS 





to preserve the tensorial character of the derived equation . 

Thus, after some manipulations we obtain the path deviation equations for an object as 
explained by MOND paradigm. 


DS^ 


1 1 

m m DS 


(38) 
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3 Path and Path Deviation of Bi-metric Theories 


3.1 Path Equation and Path Deviation of Rosen’s Approach 

Equations of path and path deviation for the Rosen version of bi-metric theory of gravity 
is derived due to the following Lagrangian: 


(39) 

where = ^ + 

By taking the variation with respect to 'I/" we obtain: 

f)T 

^ = (9.S - ( 40 ) 


dl 

dSdi>- l,.a,p)U U +{g^, ds ' 

Substituting from the above relations into the Euler-Lagrange equation; 

d dL dL 
dSd'^^ ~ ’ 

and taking into account the following conditions: 

0 , 

7/i!4|p ~ 0) 

where | is a covariant derivative with respect to . 

dU’^ 

ig,.-j,.)i^ + A;^U-un = 0, 

Multiplying equation (46) by g^^ and regarding that g’^^'y^p = 0 we obtain 


(41) 

(42) 


(43) 


(44) 

(45) 


(46) 


+ AJI.f/'C/" = 0, (47) 

which is the same equation obtained by Rosen (1940) [4], Accordingly, its corresponding 
deviation equation can be obtained by following the same technique of Bazanski approach 
by taking the variation with respect to U°‘ on (36) to obtain them - with taking into 
account the following condition to preserve its tensorial character. 


VR“ 


0 



we get after some rearrangements the following geodesic deviation eqnation: 


V52 


( 48 ) 


where is the curvatnre tensor obtained by the affine connection 7^5 [3]. 

Dne to Rosen’s approach the curvatnre tensor = 0 which reduces equation (48) 
to become 

Also, for charged objects in bi-metric theory of gravity, Falik and Rosen [33] obtained 
their corresponding held equations , which led us to introduces the following Lagrangian 
to obtain the corresponding path and path deviation equation : 


to give 


L = {g^u - 'li.uW' 


VU' 


VA 




m 


= 

VS m 

And its corresponding deviation equation becomes: 




e e 

+ -FZ^ + -(F“, - 

m w b m 


(49) 

(50) 

(51) 


Moreover, Avakian et ah [16] studied the held equations of a spinning body in the 
presence of bi-metric theory. Accordingly, we can apply the same procedure as mentioned 
in (2.1) on the following Lagrangian: 


L = (J«/J - 7«4)fl“^ + - P^ina)U‘<i‘‘S''- 


(52) 


to obtain the path equation for a spinning object in the presence of bi-metric theory : 


dU^ 


1 




dS ' 

and its corresponding spinning deviation equation: 


fiiyp/ 


(53) 


VF 2 


= (R°‘ - -\ -fi?" 

^ .pvp)^ Ld -r 2fri^ 


PU)S‘ 


up 


Vd'"' 


V5> 


Thus, if we take into consideration that P^^^ = 0, the path equation becomes: 

dU°‘ 1 

__ + = —R%pS^^U^U^ 

and its corresponding deviation equation becomes: 


(54) 


(55) 


V^dt“ 1 1 1 


VF 2 


■.pup'- 


2m 


Ds 2m 


2m 


(56) 
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3.2 Path and Path Deviation Equations of Moffat’s Approach 

Moffat [6] presented the framework of VSL satisfying bimetric theory and its causality 
to reveal the problem of dark energy due to VSL by introducing such a metric in the 
following way. 

gtiu = gtiu + Bd^(t)dy(t) (57) 

where dehnes a specihc matter metric tensor of a given matter held, B is an arbitrary 
constant has a dimension of [length]^ and chosen to be positive and ^ is a bi-scalar held 
. The inverse metrics g'^^ g^^ satisfy 

g^^g,p = s; (58) 

g^^,^ = gtiu + Bd^(l)d^(i) 

rg,p = 5 ; (59) 

Yet, the modihcation processes to control the casual propagation of the bi-scalar held led 
to redehne (57) to become: 

= g^^ + + KB^,, (60) 

where K is an arbitrary constant and is a given energy-momentum tensor to control 
the causal propagation of the biscalar held [9]. 

Consequently, we suggest the corresponding geodesic and geodesic equation owing to 
Mohat’s description of bimetric theory of gravity to be obtained by taking the variation 
with respect to on the following Lagrangian 


L = (61) 

to obtain its geodesic equation 

dU’' 

= 0, (62) 

where 

~ 2 ^ igp<^,fj- gp'^,p ~ gpp,<T)- (63) 

While, taking the variation with respect to V" , providing that = 0 to obtain its 

J-Jij 

corresponding geodesic deviation equation: 


where 




^/3-f5 




-QO. _ ■po: pQ: I ph' po; ph' pa 

■^/ 37(5 ^ ^ ^ ^7 ^ 
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3.3 Path and Path Deviation Equations of BIMOND Type The¬ 
ories 


In this section, we present the corresponding path and path deviation equation for test 
particles or spinning objects in the presence of BIMOND theories. Accordingly , it is worth 
mentioning at the beginning the above corresponding paths and there deviation equation 
in MOND paradigm to be extended in case of its BIMOND version its corresponding path 
equation becomes: 

L = + {(pa — 

where: where (p is an associated potential related to gravitational potential and 






In case of BIMOND , Milgram [14] introduced the relationship between the two affine 
connections as dehned by gfj,y and 7 ^,^ to become: 


\(X _ -pQ; pQ; 

- -L /3p - -L /3p, 


such that, 


9py,p 


QSv^pp A 9Sp^yp 


and 

~ ~ ISp^yp- 

By taking the variation with respect to 4/“ we obtain its corresponding path equation 

rlTTa 1 

While if we take the variation with respect to 17“ on the above Lagrangian we its corre¬ 
sponding deviation equation with taking into considerations the previous steps of sec( 2 . 1 ) 
to preserve its tensorial character. 

^ 

Also, in case of spinning object we introduce the following Lagrangian 

L = 9p.^^ + ^{(pa - kw + 

If we take the variation with respect to 4'“ on the above Lagrangian, its corresponding 
spinning equation becomes: 


Vf/“ 


—{R" 

2 m ^ ■ 


pup 




m 



(64) 
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and its spin deviation eqnation can be obtained by taking the variation with respect to 
t/" , with taking into consideration the steps of (2.1) to preserve its tensorial character 
we get: 


1 

= +—{R[ 




fiiyp 


1 V'h'" 

_ cup_^ _ 

2m VS 




f I i 


(65) 


3.4 Path Equations and Path Deviation of Bi-gravity Type The¬ 
ories 


Recently, Arkami et al[34] have suggested two independent metrics to explain bi-gravity 
phenomenon, 

ds'^ = g^udx^dx", 


and 

dr'^ = h^i^dx'^dx" 


Thus, in order to obtain geodesic-like equations of bi-gravity theory , we following La- 
grangian [36] : 

d dL c)L fdr 2, d dL c)L 

to give the same results as mentioned by Arkani et al (2014) 


9pv 


DU>^ 


h 


. dr , DV^ 


flU \ 


dS' Dr 


= 0 


( 66 ) 


where is an associate unit tangent vector with respect to the parameter r and 

$ is its corresponding deviation vector. 

Applying the same technique of the Bazanski approach, we obtain its deviation equations: 
to obtain: 


n2\T/« J 2 n2,T>a 

+ (—) = 0 


^dS^ 


Dt^ 


(67) 


If one considers ^ 7 ^ 0, the two metrics can be related to each other by means of a 
quasi-metric one [31]. 


9pu 9pv -|- UgUi,^ (Xhi^hgi, h)(h7)- 


( 68 ) 


Such an assumption may give rise to dehne its related Lagrangian of Bazanski’s flavor 
to describe the geodesic and geodesic deviation equation due to this version of bi-gravity 
theory. 


L = ga^U^ 




(69) 
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^/3(t 2^9 (^9<jS,I3 9S/3,(t 9l3a,s') 

and its corresponding Lagrangian: 


L 


9,uU^{ 


d'^'' 




( 70 ) 


Thus, equation of its path equation can be obtained by taking the variation respect to 
to obtain: 


DU^ 

DS^ 


(71) 


while taking the variation with respect to , with following the same steps as 
mentioned in (2.1) , we obtain its corresponding path deviation equation: 


DS^ 


■no^ 




(72) 


where 

T^cx _ ■po: _ po: _|_ per po: _ per po; 

^ ^ M^?P ' ^ ^ PP^ ^P 

Consequently, the following path and path deviation of charged and spinning objects of 
constant spinning tensor are explained as follows 


D'^9 p 

L = 

^ DS m ^ 


to give 


DU- e 


= —FPU’' 


DS m 

and its corresponding deviation equation becomes: 


~ e D^'" 

Dg2 £)5 ■■'IIP ' 


(73) 

(74) 


(75) 


where || represents the covariant derivative with respect to affine connection T^^. Also, 
the generalized path and path deviation equations for spinning objects are obtained from 
the following Lagrangian: 


Dm°‘ 1 ~ 

L = ~g,u^ + —RauuoS'^PUP^^ 


"" DS 2m 


(76) 


By taking the variation with respect T" the to obtain its corresponding path equation: 

(77) 


DTD 1 


DS 2m" 

and taking the variation with respect [/“ to obtain its path deviation equation: 
r) 2 \i>“ 1 1 D'^P 


DS^ 


DS 


(78) 
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3.5 Generalized Path and Path Deviation Equations of Bi-metric 
Theories 

Hossenfelder [35] introduced an alternative version of bi-metric theory, having two different 
metrics g and h of Lorentzian signature on a manifold M one is dehned in tangential space 
TM and the other is in its co-tangential space T*M respectively. These can be regarded as 
two sorts of matter and twin matter, existing individually , each of them has its own held 
equations as dehned within Riemannian geometry. In this part we are going to present a 
generalized form which can be present diherent types of path and path deviation which 
can be explained for any bi-metric theory which has two diherent metrics and curvatures 
as dehned by Riemannian geometry . Their Corresponding Lagrangian can be expressed 
in the following way: 

L = (79) 

By considering ^ = 0 . This will lead to two separate sets of path equations owing 
to each parameter by applying the following Bazanski-like Lagrangian: 




and 


DS 

Dt 


= 0 , 


= 0 


and their corresponding path deviation equations: 


and 


DS^ 

Dt^ 




(80) 

(81) 

(82) 

(83) 


Thus we suggest, the corresponding lagrangian to describe two independent sets of a 
generalized path and path deviation equations: 


L = 


(84) 


where, 

f, = -{eF^, + 

m 2 

and 

m 2 

By taking the variation of and respectively. 
Consequently, we obtain a set of path equations 


DU^ 

DS 


r, 


(85) 
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and 


DV^ 

Dt 


( 86 ) 


= / 


a 


and taking the variation with respect to U°‘ and using the same procedure of sec(2.1) 
to preserve their tensor character. Thus, we obtain the set of their corresponding path 
deviation equations: 


+ /“>!-'■ + g'“’hu,— (87) 

and 

a i - - D'^'" 

( 88 ) 

4 Discussion and Concluding RenicU’ks 

In this study, we have obtained the corresponding equations of path and path deviation 
equations for test particles, charged and spinning objects -constant spinning tensor- in 
different versions of Bimetric theories of gravity using a modihed Bazanski Lagrangian. 
This type of study has imposed us to determine prior to its procedure some relevant path 
and path deviation for different path equations in Weyl geometry and MOND paradigm 
to be counted as an introductory step to visualize the different stages of path and path 
deviation equations that must be included before dealing with different bi-metric theories 
of gravity. The study may give rise to search of a possible geometry able to express bi¬ 
metric theory of gravity. It can be sought that Finslerian geometry is a good candidate to 
express bi-metric theory of gravity as an extension of a Riemannian geometry [36]. In the 
mean time, path and path deviation equations using the Bazanski Lagrangian in Finsler 
geometry are in preparation [37]. Also, the above treatment of utilizing a symmetric affine 
connection, can be extended into another version of bi-metric theory of gravity following 
Einstein-Cartan geometry,as an extended approach of Drummond [38] which gives rise to 
different types of torsion and how does it propagate with respect to metric propagation 
due bimetric formalism. Finally, this work will enable us to examine, the stability of 
objects orbiting very strong gravitational field by solving the spin and spin deviation 
equations. 
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